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Abstract 

We consider the double scaling limit for a model of n non-intersecting squared Bessel 
processes in the confluent case: all paths start at time i = at the same positive value 
X = a, remain positive, and are conditioned to end at time t = 1 a.t x = 0. After 
appropriate rescaling, the paths fill a region in the fa;-plane as n — >■ oo that intersects the 
hard edge at a; = at a critical time t = t*. In a previous paper, the scaling limits for 
the positions of the paths at time t ^ t* were shown to be the usual scaling limits from 
random matrix theory. Here, we describe the limit as n — >■ oo of the correlation kernel at 
^ critical time t* and in the double scaling regime. We derive an integral representation for 

I— ' the limit kernel which bears some connections with the Pearcey kernel. The analysis is 

^ based on the study of a 3 x 3 matrix valued Riemann-Hilbert problem by the Deift-Zhou 

^ steepest descent method. The main ingredient is the construction of a local parametrix 

QO ^■t the origin, out of the solutions of a particular third-order linear differential equation, 

C ■ and its matching with a global parametrix. 
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1 Introduction and main results 

1.1 Introduction 

We considered in [25] a model of n non-intersecting squared Bcsscl paths in the confluent 
case. In this model, all paths start at time i = at the same positive value x = a > and 
end at time t = 1 at x = 0. Our aim was to study the asymptotic behavior of the model as 
n — >■ GO. 

The positions of the squared Bessel paths at any given time t G (0, 1) are a determinantal 
point process with a correlation kernel that is built out of the transition probability density 
function of the squared Bessel process. In [25] we found that, after appropriate scaling, the 
paths fill out a region in the tx plane that we described explicitly. Initially, the paths stay 
away from the hard edge at x = 0. At a certain critical time t* the smallest paths come to 
the hard edge and then remain close to it, as can be seen in Figure 1. 

In [25] we also proved the local scaling limits of the correlation kernel as n — >■ oo, that are 
typical from random matrix theory. Thus we find the sine kernel in the bulk and the Airy 
kernel at the soft edges, which includes the lower boundary of the limiting domain for t < t*. 
For t > t*, we find the Bessel kernel at the hard edge 0, see [25, Theorems 2.7-2.9] . 

In this paper we consider the critical time t = t*. We describe the transition from the 
Airy kernel to the Bessel kernel by means of a new one-parameter family of limiting kernels 
that arise as limiting kernels around the critical time. This soft-to-hard edge transition is 
different from previously studied ones in [6] or [9] , but is related to the one in [7] . 

We consider the squared Bessel process with parameter a > — 1, with transition proba- 
bality density pf given by, see [8, 23, 26], 
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Figure 1: Numerical simulation of 50 rescaled non-intersecting squared Bessel paths with 
a = 1. Bold lines are the boundaries of the domain filled out by the paths as their number 
increases. 

where la denotes the modified Bessel function of the first kind of order a, 

A remarkable theorem of Karlin and McGregor [21] describes the distribution of n inde- 
pendent non-intersecting copies of a one-dimensional diffusion process at any given time t in 
terms of its transition probabilities. In the case of the squared Bessel process, with all starting 
points at time in a > and all ending points at a later time T > in 0, the theorem implies 
that the positions of the paths at time t G (0, T) have the joint probability density 

■P(xi, . . . , x„) = ^ det[/, (j;fc)]j,fc=i,...,n <^et[gj{xk)]j,k=i,...,n (1-3) 

on (M"'")", with functions 

/2,--i(x) =x^'-ipf(a,x), j = l,...,ni := rn/21, (1.4) 

/2j(x) = x^"^p^+^(a,x), j = l,...,n2 :=n-ni, (1.5) 

gj{x) = x^-^e-^), j = l,...,n, (1.6) 

see [25, Proposition 2.1]. The constant Z„ is a normalizing constant which is taken so that 
(1.3) defines a probability density function on (M"*")". 

Formula (1.3) is characteristic of a biorthogonal ensemble [3]. It is known that (1.3) defines 
a determinantal point with correlation kernel Kn 

n 

kn{x,y) = Kn{x,y;t,T) = J2 fj(^) i^''] k,j 9k{y) (1.7) 

j,k=i 
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where [A ^]^^. is the {k,j)th entry of the inverse of the matrix 

" POO 

^= / fj{x)gk{x)dx 

-Jo J j,fc=l,...,n 



This means that 



V{xi,...,Xn) = -7 det 



{Xi , Xj ) 



i,j=l,...,n 



and for each m = 1, . . . , n — 1, 
nl /-oo /-oo 

TT / ■•■ / P{xi,...,Xn)dXm+l---dXn = det Kn{Xi,Xj) 

- my. Jo Jo I 



(n 



Ji,j=l,...,m 



(1.8) 



(1.9) 



Determinantal processes arise naturally in probability theory, see e.g. [20, 27]. The connec- 
tion with models of non-intersecting paths is well-known sec [18, Chapter 10] and references 
therein. Non-intersecting squared Bessel paths and related continuous models with a wall 
are studied in [22, 23, 24, 26, 28]. Non-intersecting discrete random walks with a wall are 
considered in the recent papers [4, 5, 7, 29]. 
As in [25] we introduce a time rescaling 



t ^ 



2n 



1 

2^ 



and we consider the rescaled kernels 



Kn{x,y;t) = e-^^--yy^'-'^Kn(^x,y 



t 1 
' 2^' 2^ 



x,y>0, <t <1, 



(1.10) 



that depend on the variable t. The prefactor e "'^^ y)/i^ *) does not affect the correlation 
functions (1.9). We define wi^n, W2,n on [0, oo) by 



^2,n(x) = a;("+i)/2 



2n\/ax 



exp 



nx 



t{l-t) 



t 



l-a+l 



2n\fax 



t 



(1.11) 



as in [25, equation (2.20)]. 

Then the kernel (1.10) is expressed in terms of a RH problem. Indeed we have 



Kn{x,y;t) = 



1 



2iTi{x — y) 



(0 «;i,n(y) w2,n{y))y+\y)y+{^) 



(1.12) 



where Y is the solution of the following matrix valued Riemann-Hilbert problem, see [25]: 
RH problem 1.1. Find y : C \ M ^ C^^^ ^^^^ ^^^^ 

1. Y is analytic in C \ [0, oo). 

2. On the positive real axis, Y possesses continuous boundary values i+ (from the upper 
half plane) and YL (from the lower half plane) , and 



'I Wi^n{x) W2,n{x) 

F+(x) = F_(x) ( 1 
,0 1 



X > 0, 



(1.13) 
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3. Y{z) has the following behavior at infinity: 

Y{z)= [l + 0{-^^ \ |, z^oo, zeC\W, (1.14) 










— ni 











where ni = [n/2] and n2 = [n/2\. 
4. Y{z) has the following behavior near the origin, as 2; — >■ 0, 2; G C \ [0, 00), 



'1 h{z) 1 
Y{z) = \1 h{z) 1 I , with h{z) = < 
,1 h{z) 1 



f \zf, if -1< a < 0, 

log|z|, if a = 0, (1.15) 
1, if < a. 



The O condition in (1.15) is to be taken entry-wise. 

This RH problem has a unique solution given in terms of multiple orthogonal polynomials 
for the modified Bessel weights (1.11). 

It was proven in [25, Proposition 2.3 and Theorem 2.4] that in this scaling there is a 
critical time 

f = (1.16) 

depending only on the starting value a. For every t G (0, 1), we have that 

lim —Kn(x,x;t) = p(x;t) 

n->-oo n 

exists, where the limiting density p{x;t) is supported on an interval \p{t),q{t)] with p(t) > 
ii t < t* and p{t) = if t > t*. The results of [25] were obtained from a steepest descent 
analysis of the above RH problem for values oit^t*. In this paper we develop the steepest 
descent analysis at the critical time. 

1.2 Statement of results 

The main result of our paper is the following theorem. 

Theorem 1.2. Let Kn he the correlation kernel (1-12) for the positions of the rescaled non- 
intersecting squared Bessel paths starting at a > at time and ending at zero at time 1. 
Let t* = a/{a + 1) as in (1.16) and 

Then we have, for every fixed r G M, and x,y > 0, 

lim 4j:-Kj^,^;t*-^) = K^"\x,y;T), (1.18) 



where K^"^ is the kernel 



KZ^'ix, y- T) = [ [ *le-'/t+l/i2t')-r/s-l/i2s-)^xt-ys dtd^ (^ ^g) 

[zTn) 7ter Jse's s — t 

The contours T and S in (1.19) are as in Figure 2. The fractional powers s" and in (1.19) 
are defined with a branch cut on the positive semi-axis, i.e., < args,argi < 2%. 
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Figure 2: The contours of integration F and S used in the definition of the critical kernel 
(1.19). The contour F consists of a closed loop in the left half-plane tangent to the origin and 
is oriented clockwise. The contour S is an unbounded loop oriented counterclockwise and 
encircling F. 

We prove Theorem 1.2 by an asymptotic analysis of the RH problem 1.1 by means of the 
steepest descent analysis of Deift and Zhou, as we did in [25] for the non-critical times. 

At a certain stage in the analysis we have to construct a local parametrix at the origin 
X = (the hard edge). This was done in [25] with the Besscl parametrix. We also had to 
construct an Airy parametrix at another position (a soft edge). In the critical case that we 
are considering in this paper this other position coincides with the origin. The coalescing of 
the soft edge with the hard edge leads to the construction of a new local parametrix at the 
origin. The construction uses a new model Riemann-Hilbert problem that we describe in the 
next subsection. The functions that appear in the model RH problem ultimately lead to the 
expression (1.19) for the limiting kernels. 

1.3 Riemann-Hilbert problem 

The model RH problem is defined on the contour S$ shown in Figure 3. It consists of the six 
rays argz = 0, ±7r/4, ib37r/4, oriented from left to right. 

RH problem 1.3. Let a > — 1 and r G C. The RH problem is to find $q = $Q(-;r) : 
C \ S$ ^ C3^3 s^gj^ that 

1. $a is analytic in C \ S$. 

2. $a has boundary values on each part of S$ \ {0} satisfying 



(1.20) 



where the jump matrices are shown in Figure 3. 

3. Let CO = e^'^*/^ ^nd 



A; = 1,2,3. 



(1.21) 
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Then as ^ — ^ cx), we have 



$^(z;r) = ^^^ I 1 I I 1 1 l| I e-"^*/3 Q 



-a/3 



Z 











1 


i)( 




a; 





1; 



(/ + 0(z-^/^)) e^2(2;r) , lmz>0, (1.22) 
^ ^ ^ e^3(.;r) ' 



and 

-a/3 



i/''^ \ _^ i\ /g-aTTi/s Q 



$^(^;r) = ^^^ I 1 I ( 1 -1 1 I I e"'^'/''^ 



V3 



Z"^/V 1/ \ ly 

/e^2(^;r) \ 

(/ + C»(z-i/3)) e^i(^;^) , lmz<0, (1.23) 

^ Mo e^3(.;r) / 



4. As z — )• we have 



Iz'^ 0> 

$a(z;r) 0|=O(l), < |argz| < 7r/4, (1.24) 

Vo \j 

(I 0^ 

$a(-2;T) z° 0|=O(l), 7r/4 < larg^l < 37r/4, (1.25) 

Vo 1. 

/I 

$a(z;r) I =0(1), 37r/4 < |argz| < TT. (1.26) 

VO 

Note that the parameter r appears in (1.21) and in the asymptotic conditions (1.22) and 
(1.23) of the RH problem. 
We prove the following. 
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Theorem 1.4. Let a > —1 and r G C. The RH problem 1.3 for has a unique solution 
with 

det ^ a{z; t) = z-"", zeC\T.^. (1.27) 
The critical kernel (1.19) satisfies 



i^r*(^,y;r) = ^^.J_^^ (-1 1 0)$-V(2/;r)$.,+ (x;r)|l|, (1.28) 



for x,y > and r G 



The uniqueness statement in Theorem 1.4 follows from standard arguments where one 
first proves (1.27). The existence of a solution follows from an explicit construction of ^'a, 
given in Proposition 5.2, in terms of solutions of the third order ODE 

zp'" + ap" -Tp' -p = (}. (1.29) 

A particular solution of this equation is given by 

where F is the closed contour shown in Figure 2. 

The inverse matrix is built out of solutions of the adjoint equation 

z(^" + (3 - a)g" - rg' + g = (1.31) 

which has the special solution 



where S is also shown in Figure 2. 

In terms of these functions the kernel (1.19), (1.28), can also be written as 

2ni{x-y)K^"\x,y;T) 

= [q"{y) -{a- 2)q'{y) - Tq{y)] p{x) + [-yq'{y) + (a - l)q{y)] p'{x) 

+ yq{y)p"{x). (1.33) 

For y = X the right-hand side of (1.33) is the bilinear concomitant [1, 19] which is constant 
for any two solutions p and q of the differential equations (1.29), (1.31), and which turns out 
to be zero for the two particular solutions (1.30) and (1.32). 

Remark 1.5. There are solutions of the differential equations (1.29) and (1.31) that can be 
written as integrals of Bessel functions. In particular, we have that 



/■+00 

p{z) = z-"/^ u"/2e-™-«'/2j^(2V^)d« and 

^0 

/+ioo 
^-a/2^rv+v /^J^(2^)dv, 
-ioo 



(1.34) 
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are solutions of (1.29) and (1.31), respectively, where Jq is the Bessel function of the first 
kind of order a. 

Based on a similarity with formulas by Desrosicrs and Forrester [17, Proposition 5] for a 
perturbed chiral GUE, we suspect that it should be possible to write an alternative expression 
for the critical kernel in (1.19) in terms of the functions (1.34), namely 

dudv 



f f ™j„(2Vxu)J„(2Vy^) 



u — V 



Unfortunately, we have not been able to make this identification. Observe however that 
for a = —1/2 the double integral above reduces to the so-called symmetric Pearcey kernel 
^(ci; 172; r/), with ai = x^/\/2, 02 = y^/\/2, 77 = \/2r. The correlation kernel 

^{av,a2;v) = ^. I I 6"''-'+''-'--'+"' cos(aix) cos(a2u)^^^, 

* Ju&C Jx€lS.-i- U — X 



where C is the contour in C consisting of rays from ooe*'^/^ to to ooe"*'^/^, was introduced 
by Borodin and Kuan in [7]; the authors point out the possible connection with the non- 
intersecting Bessel paths in the critical regime, as it seems to be the case. 



2 First and second transformation 

The steepest descent analysis consists of a sequence of transformations 

Y ^ X S ^ R 



which leads to a RH problem for R, normalized at infinity and with jump matrices that are 
close to the identity matrix if n is large. 

We start from the RH problem 1.1 for Y, stated in the introduction. The RH problem 
depends on the parameters n and t. We assume that n is large, and t is close to the critical 
value t*. Eventually we will take the double scaling limit 

n ^ 00, t ^ t*, such that \fn{t*' —t) = c*t remains fixed. (2.1) 

But throughout the transformations in Sections 2-7, we assume that n and t are finite and 
fixed. 

The first transformation is the same as in [25] . 



2.1 The first transformation 

The first transformation Y X is based on special properties of the modified Bessel func- 
tions that appear in the jump matrix (1.13) via the two weights (1.11). The result of the first 
transformation will be that the jump matrix on [0, 00) is simplified at the expense of intro- 
ducing jumps on (—00,0) and on two unbounded contours that are shown in Figure 4. 
The contours are the boundaries of an unbounded Icnse around the negative real axis. 

Here and in the sequel, Eij denotes the 3x3 elementary matrix whose entries are all 0, 
except for the (i,j)-th entry, which is 1. 
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Figure 4: Contour Sx 



. U and jump matrices Jx in the RH problem for X. 



Definition 2.1. We let yi{z) = z("+i)/2l^+i(2/^) and 2/2(2) = z("+i)/2iC„+i(2^) where 
ii'a+i is the modified Bessel function of second kind of order a + 1. Then we define X in 
terms of Y as follows 




-z "yi 



\ (I 





1^ 

















-2TTi 



X < 



1 — e Z i^23 



I + e-^'^'z"" 



for z in the upper part of the lens, 
£'23; for z in the lower part of the lens, 
elsewhere. 



(2.2) 



where Ci is some constant matrix, see [25, Equation (3.12)] for its definition. 

Then X is the unique solution of the following RH problem, see [25, Section 3] for details. 
RH problem 2.2. 

1. X is defined and analytic in C \ 'Ex where Sx 

2. On Sx we have the jump 



.UA±. 



X+ = X^Jx 
where the jump matrices Jx are as in Figure 4. 
3. As z — )• 00 we have 



(2.3) 



X{z) 



I + O 











^{-l)"/4 













a/2 







^-a/2 




/l 

71 71^ 

\' 75^ 75, 


Q Q ^—n 1 2 ^n^faz 1 1 ^ 



(2.4) 
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4. X{z) has the same behavior as Y{z) near the origin, see (1.15), as z ^ from outside 
the lens around (— oo,0]. If 2; ^ within the lens around (— oo,0], then 



X{z) = { 



(1 1' 

o 1 1 

\l \z\- I, 

(1 logl^l log|z| 
1 logl^l logl^l 
1 log|z| log 1^1 
/I 1 \z\ 
C I 1 \z\ 
Vl 1 \z\ 



if a < 0, 
if a = 0, 
if a > 0. 



(2.5) 



2.2 The Riemann surface 

In the second transformation we are going to use certain functions that come from a Riemann 
surface. In [25, Section 4] we used the Riemann surface associated with the algebraic equation 



z = 



1 - K 



k = {l-t){t- a{l-t)). 



ai-t{i-t)cr 

This equation was derived from a formal WKB analysis of the differential equation 

2nz 



(2.6) 



zy'"{z)+[{2 + a) 



+ 



n^z n{n — a — 2) an^ 



t2{i _ tf + i(l - t) 



n 



^2 jy(^) i2(i-f)2 



y(z) = 0, (2.7) 



see [11] and [25, Equation (2.21)], that is satisfied by the multiple orthogonal polynomials 
associated with the weights (1.11). 

There are three inverse functions to (2.6), which behave as z — >■ 00 as 



Cl(2) 



2 / ' 



C2iz) = 



C3(^) 



t{l - 1) tz^^ 2z 



1 ^ 



+ 



^3/2 
1 



(2.8) 



t{\-t) tz^^ 2z V^^/^ 
At critical time t = t*, we have = and equation (2.6) reduces to 

1 



C(1-C*C)2' 



C* =t*{\-t*). 



(2.9) 



Then, the corresponding Riemann surface has two real branch points, and q* = 27c* /4 > 0, 
being degenerate (of order 2), and q* being simple. The point at infinity is also a simple 
branch point of the Riemann surface. 
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l/3c 



1/c 



w 



Figure 5: Plot of function z = z{w) given by (2.10) for G M. 



In the present paper, we want to work with a Riemann surface TZ with a double branch 
point, even if t 7^ t*. Following the approach of [2] we do not consider (2.6) if t ^ t* but 
instead consider a modified equation 

z=—^-—^, (2.10) 

with c some positive number and w a new auxiliary variable. 
The Riemann surface has simple branch points at 



and at infinity, and it has a double branch point at 0. The sheet structure of TZ can be readily 
visualized from Figure 5 and is shown in Figure 6. 

The sheets TZi and 7^2 are glued together along the cut Ai = [0, q\ and the sheets 7^3 and 
7^2 are glued together along the cut A2 = (0, 00] in the usual crosswise manner. There are 
three inverse functions tyj^, A; = 1, 2, 3, that behave as 2; — >■ 00 as 



and which are defined and analytic on C \ Ai, C \ (Ai U A2) and C \ A2 respectively. 

The algebraic function w = w{z) in (2.10) gives a bijection between the Riemann surface 
71 and the extended complex tu-plane. Figure 7 represents this mapping, along with the 



q = 27c/4 



(2.11) 




(2.12) 
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Figure 6: The Riemann surface TZ used in the steepest descent analysis. The origin is a double 
branch point for every t. 

domains itj = Wj{Tlj), j = 1,2,3 (the images of the corresponding sheets of TZ) and the 
points 

Wq = W2{q) = — > 0, Woo= W2(po) = - > 0. 

ic c 
We point out that W2+{^j), j = 1,2, are analytic arcs that extend to infinity in the upper 
half plane, while W2-{Aj), j = 1, 2, are in the lower half plane. 

2.3 Modified ( functions 

We next define modified ^-functions with the same asymptotic behavior as z — >■ oo as given 
in (2.8) up to order ©(z"^/^). 

Definition 2.3. For /c = 1, 2, 3 we define 

Ck = Wk+pwl, (2.13) 

where Wk is given by (2.12) and 

(i-t)2/va /« 2rY , 

Note that for t = t* , the critical time, we have c = c* = t*{l — t*) and p = 0, so that we 
recover in this case the equation (2.9) and the (^-functions defined in (2.8). 

Lemma 2.4. For c andp given by (2.14), the asymptotic behavior of functions (k> k = 1,2, 3, 
defined in (2.13), as z ^ oo, is given by (2.8). 

Proof. This follows from direct computations using (2.12) and (2.13). □ 

In what follows we need the behavior of the ^-functions (2.13) near the origin. The 
following lemmas are analogous to Lemmas 3.2 and 3.3 in [2]. We put before. 
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Lemma 2.5. There exist analytic functions fi and gi defined in a neighborhood Ui of the 
origin such that with z £ Ui and k = 1,2,3, 

^ (^) = W'^z-^'^h iz) + uj^z^^gi {z) + I for Im z > 0, 

\u''z-^l'^fi{z)+uj'^^z^l^gi{z) + ^^ /or lmz<0, ^ " ' 

In addition, we have /i(0) = c~^/^, and fi{z) and gi{z) are real for real z G Ui. 

Proof. We put z = and w = (xy)"^ in (2.10) and obtain 

ixy-cf-y' = 0. (2.16) 

It has a solution y = y{x) that is analytic in a neighborhood Ui of and y(0) = c^^^ > 0, 
y'(0) = — |c^/^ < 0. Then, we can write 

xw{x) = l/y{x) = fi{x^) + x^gi{x^) + xhi{x^), (2.17) 

with /i, gi and /ii analytic in Ui and /i(0) = c^^/'^ > 0. Plugging this into (2.10), we find 
after some calculations that 

fiiz)giiz) = ^, c^ff(z)-l + ^^z + c^gfiz)z^ = 0, (2.18) 

and hi{z) = 2/(3c). Hence, there is a solution w = w{z) of (2.10) with 

w{z) = z-'/^fi{z)+z^/^gi{z) + for zGUi\ (-oo,0], 

3c 

where we take the principal branches for the fractional powers. This solution is real for z real 
and positive, thus it coincides with w^lz), which proves (2.15) for k = 3. Considering the two 
others solutions of (2.16), we get the expressions (2.15) for /c = 1, 2 by analytic continuation. 
Since y{x) is real for real x, (2.17) implies that fi{z) and gi{z) are also real when z is real. □ 

The next lemma describes the behavior of the (^-functions at the origin. 

Lemma 2.6. There exist analytic functions f2 and g2 defined in a neighborhood U2 of the 
origin such that for z £ U2 and k = 1,2, 3, 

ru.2fc^-l/3/2(z)+^'=Z-2/3^2(^) + 3^(^ /orlmz>0, 
|^^'^-l/3/2(z)+^2fc^-2/3^^(^)^_2_ forluiZ<0. ^" ' 
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In addition, we have 

/2(0) = c-2/3 + ^) , 52(0) = (2.20) 

and f2{z) and g2{z) are real for real z E 1/2- 

Proof. The proof follows from (2.13) and the previous lemma. It suffices to compute w^{z) 
by using (2.15) and the first identity in (2.18). Then, (2.19) follows from (2.13) if we set 

f2iz) = (1 + f ) /i(^) +P^9fi^l 92{z) = (1 + f ) ^51 (^) +Pf!{z), (2.21) 

and (2.20) follows from the value of /i(0) given in Lemma 2.5. The functions f2{z) and g2{z) 
are real for real z e U2 since fi{z) and gi{z) are real for real z eUi. □ 

2.4 The A-functions 

We next define the A-functions as anti-derivatives of the modified ^-functions (2.13). 
Definition 2.7. We define for A; = 1, 2, 3, 

Xk{z) = r Ck{s)ds, (2.22) 

where the path of integration starts at on the upper half-plane (which is denoted by 0_|_) 
and is contained in C \ (—00, q] for A; = 1, 2, and in C \ (—00, 0] for A; = 3. 

By construction the functions Ai and A2 are analytic in C \ (— oo,gr] while A3 is analytic 
in C \ (—00, 0]. Prom Lemma 2.4 and (2.22), it follows that, as z ^ 00, 



Xi{z) = logz + £i + 



r 

z ^ 

X / X ^ 2VazV2 1 / 1 \ 



where ^fe, k = 1,2,3, are certain constants of integration, and log z is the principal branch of 

the logarithm which is real on (0, +00). Using the structure of the Riemann surface Tt and 
the residue calculation based on the expansion of C,\ at infinity, see (2.8), we conclude that 

Ai_(0) = -27rz, A2-(0) = 27ri, A3-(0) = 0. (2.23) 

Moreover, the following jump relations hold true: 

Ai±(a;) = A2=F(a;) - 27ri, A3+(a;) = A3_(x), xGAi = (0,(j), 

Ai+(.x) = Ai_(x) + 27ri, a; G A2 = (-00, 0), (2.24) 

A2+(a;) = A3_(,x), A2-(x) = A3+(a;) + 27ri, x G A2. 

The behavior of the A-functions at the origin follows from Lemma 2.6 and (2.22). 
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Re>.^ = ReXj 

Rek^ = ReXg 

ReX^ = ReK 



-4 



Figure 8: Curves ReAj = ReA^ at the critical time (here a = 1 and t = t* = 0.5). 

Lemma 2.8. There exist analytic functions /s and defined in a neighborhood Us of the 
origin such that for z £ U3 and k = 1,2,3, 

X = ilu;''z'/^f3iz)+u;'z'/^gsiz) + 3^ forlmz > 0, 

|Afe_(0) + |a;V/3/3(-^)+c^'V/353(z) + 3^ forlmz<Q, ^' ' 

with Ajt_(0) given by (2.23). In addition, we have 

fsiO) = /2(0) = (1 + ^) , 53(0) = 352(0) = 3pc-V3, (2.26) 
and both fs{z) and gs{z) are real for real z E U3. 

Proof. Integrating (2.19) and taking into account that Afe+(0) = 0, we get (2.25) and (2.26). 
The fact that fsiz) and gsiz) are real for real z G C/3 also follows from Lemma 2.6. □ 

The functions /a and 53 depend on t. We write f3{z; t) and 53(2;; t) if we want to emphasize 
the dependence on t. 

In order to control the jumps in the different RH problems that we are going to consider in 
the sequel, we need to compare the real parts of the A-functions. Figure 8 shows, at the critical 
time t = t*, the curves in the complex plane where the real parts of the A-functions are equal. 
These curves are critical trajectories of the quadratic differentials {Cj{z) — Qkiz))"^ dz"^ . The 
curve ReA2(2) = ReA3(z) consists of the negative real axis. The curve ReAi(2;) = ReA3(2;) 
consists of two trajectories emanating from the origin, symmetric with respect to the real 
axis, and going to infinity. Finally, the curve ReAi(2;) = ReA2(z) consists of Ai along with 
two symmetric trajectories emanating from the branch point q* . 
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Figure 9: Curves ReXj = Re near the origin, before the critical time (here a = 1, i = 0.3 < 
t* = 0.5). 

At a non-critical time t ^ t*, the configuration of the curves remains the same, except in 
a small neighborhood of the origin. Figures 9 and 10 show the new configurations in such a 
neighborhood. When t < t*, the function Ci{z) — C2(^) has an additional zero on Ai, which 
causes the appearance of a new loop around the origin where ReAi(z) = Jie\2{z). Also, the 
curve ReAi(2;) = ReX^lz) is shifted to the left and becomes a continuation of the loop as it 
intersects the negative real axis. Similarly, when t > t*, the function C2{z) — Csi^) l^as an 
additional zero on A2, which causes the appearance of a new loop around the origin where 
ReA2(-2) = ReA3(2;). Now, the curve ReAi(2;) = ReA3(z) is shifted to the right and becomes 
a continuation of the loop as it intersects the positive real axis. In both cases, when t tends 
to t*, the loop shrinks to the origin. 

Concerning the relative ordering of the real parts in a neighborhood of the real axis, the 
following lemma holds true. 

Lemma 2.9. (a) For z G (g, +00), we have ReAi(z) < ReA2(2;). 

(b) The open interval (0, q) has a neighborhood U\ in the complex plane such that for z G 
Ui \ (0, q) and outside of the additional loop around when t ^ t* , we have Re A2(2:) < 
ReAi(z). 

(c) The open interval (—00, 0) has a neighborhood U2 in the complex plane such that for z G 
C/2\(— 00, 0) and outside of the additional loop around when t ^ t* , we have Re A2(^) < 
Re A3 (2) . The neighborhood U2 is unbounded and contains a full neighborhood of infinity. 

Proof This is similar to the proof of [25, Lemma 4.3]. When t ^ t* , only the ordering of the 
real parts inside the additional loop is modified. See also [2, Lemma 4.2]. □ 
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0.01 
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ReX^ = Re 

Rek^ = Re 
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-0.025 



Figure 10: Curves Re Xj = Re near the origin, after the critical time (here a = l,t = 0.7 > 
t* = 0.5). 

In the double scaling limit (2.1) that we are going to consider, we have that t — t* is of 
order n~^/^ as n — > oo. Then, the special ordering of the real parts of the A-functions inside 
the loop will not cause a problem because the (shrinking) disk around the origin where we 
are going to construct the local parametrix will be big enough to contain the loop for n large, 
see the proof of Lemma 6.3 below. 



2.5 Second transformation of the RH problem 

The goal of the second transformation X t-^ U is to normalize the RH problem at infinity 
using the functions Xj from Section 2.4. 



Definition 2.10. Given X as in (2.2) we define 



/g-nAi(z) 



Uiz) = C2X{z) 



\ 











-n{X2{z)~ 



t(i-t)' 








e-'»(^3(^)-^) J 

where C2 is some explicit constant matrix, see [25, Equation (4.19)]. 



(2.27) 



As a consequence of the assertion (c) in Lemma 2.9 we may (and do) assume that the 
contours A^, defined above (and depicted in Figure 4) lie in the neighborhood U2 of A2 where 
Re (A2 — A3) < (except when it intersects the small loop near when t ^ t*, see Figures 9 
and 10). 

Using the jump relations (2.24) and other properties of the A functions we find that U 
solves the following RH problem, see [25, Section 4] for details. 
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RH problem 2.11. The matrix-valued function U{z) defined by (2.27) is the unique solution 
of the following RH problem. 

1. U{z) is analytic in C \ Sf/ where Sc/ = Sx = K U A^. 



2. On T,if we have the jump 

U+{z) = U.iz)Juiz), 
with jump matrices Ju{z) given by 

/I 

Ju{x) = -|a;|- 
\0 

gn(A2-Ai)_(x) Q 

L 

Juix) = (/ + x"e'*(^i-^2)(^)£;i2) , 

3. As z ^ oo we have 



zeT, 



u 



X G A2 = (-00,0), 

XG Ai = {0,q), 
X G {q,oo), 



G A^. 



(2.28) 

(2.29) 

(2.30) 

(2.31) 
(2.32) 



U{z) 



I + O 










1/4 








^-1/4 



/I 





\ 





1 

V2 V2' 



1 

V2 . 







-a/2 



(2.33) 



4. U{z) has the same behavior as X{z) at the origin, see (1.15) and (2.5). 



It follows from Lemma 2.9 that in the double scaling limit (2.1) the jump matrices in 
(2.31) and (2.32) tend to the identity matrix / as n — )■ 00 at an exponential rate, except when 
z lies inside the small loop around the origin, see Figures 9 and 10. Moreover, (A2 — Ai)_)_ = 
— (A2 — Ai)_ — A-Ki is purely imaginary on Ai, so that the first two diagonal elements of the 
jump matrices in (2.30) are highly oscillatory if n is large. 

3 Third transformation of the RH problem 

The third transformation U ^ T consists in opening a lens around Ai, see Figure 11. It 
transforms the oscillatory entries in the jump matrix on Ai into exponentially small off- 
diagonal terms. Following [25, Section 5], we define T as follows. 

Definition 3.1. We define 

T{z) = U{z) (l T z-°e"(^2-^i)(^)E2i) , (3.1) 

for z in the domain bounded by A^ and Ai (the shaded region in Figure 11), and we define 

T{z) = U{z) (3.2) 

for z outside of the lens around Ai. 
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Figure 11: Opening of lens around Ai. The contour Sj- = M U U is the jump contour 
in the RH problem for T 



Let Er = M U Af U Af be the union of the contours depicted in Figure 11. Then, 
straightforward calculations yield that the matrix-valued function T is the solution of the 
following RH problem. 

RH problem 3.2. The matrix-valued function T(z) satisfies 

1. T is analytic in C \ S^. 

2. On Sj' we have 



T-Jt 



where 



Jt{x 

Jt{z 

Jt{x 

Jt{z 
Jt{x 




\x\ 




J _|_ ^±am ^-a ^n{\2~\-i){z) 



E23, 




X G A2, 



X G Ai, 

z £ Af, 

X G {q, +00). 



(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 
(3.8) 



3. As z — )• 00, we have 



T{z) 



I + O 










1/4 








^-1/4 



/I 



— —I 



4. For — 1 < Q < 0, T{z) behaves near the origin like: 



'10 

,0 



'1 1^ 
T{z) = \l 1 I , as z ^ 0. 

.1 kl" I, 



(3.9) 
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For a = 0, T{z) behaves near the origin hke: 



'1 log|z| 1^ 

T{z) = O \ 1 log \z\ 1 I , as z — >■ outside the lenses around A2 and Ai, 
,1 log b 1, 



and 



/l log|z| log|z| 



T{z) 



O 



as 2; ^ inside the lens around A2 , 



O 



as 2; ^ inside the lens around Ai. 



1 log l^l log \z 
\1 \og\z\ \og\z 
''logl^l log 1 2; I 1^ 
log|z| log 1 2; I 1 

yiogi^i log 1 2; I 1/ 

For a > 0, T{z) behaves near the origin like: 



/I 1 1^ 

T{z) = \l 1 l|, as2;^0 outside the lenses around A2 and Ai, 

VI 1 i> 

and 



T{z) = 



1 


1 


\z 


1 


1 


\z 


1 


1 


\z 


\z 


—a 


1 


\z 


—a 


1 


\z 


—a 


1 



as 2; ^ inside the lens around A2, 
as 2; ^ inside the lens around Ai. 



5. T{z) remains bounded as 2; — >■ 

4 Global parametrix 

In the next step we ignore the jumps on A^ and A^ in the RH problem for T and we consider 
the following RH problem. 

RH problem 4.1. Find AT^, : C \ (-00, q] C^^^ such that 

1. Na is analytic in C \ {—oo,q\. 

2. Na has continuous boundary values on (—00, 0) and (0, q), satisfying the following jump 
relations: 



) 







—a 








) 




















\x\ 










xe {0,q), 
X G (—00, 0). 



(4.1) 
(4.2) 
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3. As z — >■ oo, 



10 \ 

72^ 72/ 



'10 
z"'/^ 



,0 



,-a/2 ^ 



4. As z — > we have 



(4.3) 



(4.4) 



5. As z ^ we have 



/I 0^ 

2"/^A^a(z) z-'^ 

\0 1> 

where is a rank one matrix. 



(4.5) 



This RH problem is solved as in [25, Section 6] in terms of the branches Wk, k = 1,2,3, 
of the algebraic function w, defined by (2.10), (2.12). 



Proposition 4.2. (a) The solution of the RH problem 4-1 for a = is given by 

/Fi{wi{z)) Fi{w2{z)) Fi{w3{z))\ 

No{z) = F2{wi{z)) F2{W2{Z)) F2{W3{Z)) , 

\F3iw1iz)) F3{w2{z)) Fsiwsiz))/ 



(4.6) 



where 



„ / X (w — Woo)'^ „ , s w(w — W*) „ , s w(w — Woo) , . _N 

with w* ^ Woo, Ki, K2, K3 certain non-zero constants that depend on a and t, and 

D{w) = {w - Wq){w - Woo), (4.8) 

The square root D{w)^/'^ in (4.7) is defined with a cut along W2-{^i) U W2-{1^2), such 
that it is positive for real w > Woo ■ 



(b) The solution of the RH problem 4-1 for general a is given by 

No,{z) = CaNoiz 



e"'^2(2) 



(4.9) 



e"<^3(^) 



where No{z) is given by (4.6), Cq, is a constant matrix that depends on a, t and a, 
such that detCa = 1 (see [25, Equation (6.I4)] for details), and the functions Gj{z) 
are given by 

Gj{z) = rj{wj{z)), zeTZj, j = 1,2,3, (4.10) 
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with 



ri{w) = log(l — cw), w £ TZi, 

r2{w) = — logw; — log(l — cw), w £ Tl2, 
fsiw) = log(l — cw) + iir, w E TZ^. 



(4.11) 



The branches of the logarithms in (4.11) are chosen so that log(l — cw) vanishes for 
w = and has a branch cut along W2-{A.2) (cf. Figure 7), and logw is the principal 
branch with a cut along (— oo,0]. 

Proof. The fact that satisfies items 1., 2., 3., and 4. of the RH problem 4.1 is proved as 
in [25, Section 6]. 

Prom (4.6), (4.7), and the behavior of the functions Wk at as given in Lemma 2.5 we 
obtain 



A^o(^) 







C-2/3 


K2 


< 








C-2/3 









a;2 Lo 1)^-^3 + 0(1)^ z^O,Imz>0, 



z ^ 0, Imz < 0, 



(4.12) 



where, for Im z < 0, there is a minus sign in the second entry of the rowvector (a; — l) 
because of the choice for the branch of the square root D{w)^/^ used in (4.7). This proves 
(4.5) for the case a = 0. 

For general a, we use that by (4.10) and (4.11) we have that 

^G,iz) ^ c»(z-V3), e«^(^) = 0(^2/3), e«3(.) = c»(2-V3) 

as 2; ^ 0. Then using (4.9) we find (4.5) also for this case. □ 
It will be convenient in what follows to consider besides also the matrix valued function 



NJz) 




(4.13) 



which already appeared in (4.5). 

Lemma 4.3. The solution Na of the RH problem 4-1 is unique and satisfies 

det Na{z) = det Na{z) = 1, zeC\{-oo,q]. 

where Na is given by (4.13). 
As z ^ we have both 

Na{z) = O {\z\-^l^^ , and N^iz)'^ = O {\z\-^/^^ . 



(4.14) 
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Proof. Prom the jump condition in the RH problem 4.1 it follows that det Na{z) has an 
analytic continuation to C \ {0,q}. The isolated singularities are removable by (4.4) and 
(4.5), where it is important that in (4.5) is of rank one. Thus del Na{z) is an entire 
function, and since it tends to 1 as ^; ^ oo by (4.3) we conclude from Liouville's theorem 
that det Na{z) = 1. The uniqueness of the solution of the RH problem 4.1 now follows with 
similar arguments in a standard way. 

From the definition (4.13) we then also find that det Na{z) = 1. The behavior (4.14) for 
both Na and its inverse finally follows from the condition (4.5) in the RH problem, and the 
fact that Na has determinant one. □ 

We need two more results that will be used later in Section 7. 

Lemma 4.4. For Nq defined in (4.9), we have 

fl 0\ 

Nq^ = N^\0 zeC\{-oo,q], (4.15) 

\0 -i / 

where the superscript T denotes the matrix transpose. 

Proof. Observe from (4.1)-(4.2) that Nq and Nq'^ have the same jumps on Ai and A2, so 
that NqNq is analytic in C \ {0, q}. The singularies at and q are removable because of (4.4) 
and (4.5), so that NqN^ is entire. By (4.3), 

fl 0\ 

No{z)N^{z) = i +O - , z^oo, 
\0 i 0/ 

and the assertion (4.15) follows by Liouville's theorem. □ 

As a consequence, we obtain the following corollary. 
Lemma 4.5. The constants Kj from (4.7) satisfy the relation 

Kl - 2iK2K3 = 0. (4.16) 

Proof. Prom (4.15) we obtain 

fl 0\ 

N^{z) -i\ No{z) = I. 
\0 -i / 

Then insert the behavior (4.12) for both Nq{z) and No{z) and observe that the coefficient of 
^-2/3 jnust vanish. This yields 

/I \ fK,\ 
{Kl i^3) -i K2 = 0, 

\o -i 0/ yKj 

which is (4.16). □ 
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5 Local parametrices 

The next step is the construction of local parametrices around the branch points q and 0. We 
shall be brief about the local parametrix P around q. The main issue will be the construction 
of the local parametrix Q around the origin. 

5.1 Parametrix P around q 

We build P in a fixed disk D{q*,rq) around q* = 27c* /A with some (small) radius > 0. 
For t sufficiently close to t* we then have that q is in this neighborhood, and we ask that P 
should satisfy. 

1. P is analytic on D{q*,rq) \ Ey; 

2. P has the same jumps as T has on D{q*,rq) fl St, see (3.6)-(3.8); 

3. as n — >■ oo, 

P(z) = Na{z) {I + O (n"^)) uniformly for \z - q*\ = Vq. (5.1) 

The construction of P is done in a standard way by means of Airy functions, see [13, 14, 15] . 
We will not give any details. 

5.2 Parametrix Q around 0: required properties 

The construction of the parametrix at the origin is the main novel ingredient in the present 
RH analysis. A similar problem has been previously solved in [2, Section 8], which serves as 
an inspiration for the approach we follow here. 

We want to define a matrix Q in a neighborhood D(0, ro) of the origin such that 

1. Q is analytic on £)(0, ro)\ST, where has been defined in Section 3, see also Figure 11; 

2. Q has the same jumps as T has on n D(0, ro), see (3.4)-(3.7). That is, we have 

Q+{z) = Q.{z)Jq{z), ^GETnL'(0,r), (5.2) 

where Jq is given by 

/I \ 

Jq(x)= -l^r^^ , xG A2nL>(0,r), (5.3) 

\0 / 

Jq(z) = / + e±"'^V°^e"(^2-^3)(^)^23, z G A^nL>(0,r), (5.4) 

/ x° 0\ 

Jq{x)=\-x-'^ , xGAinD(0,r), (5.5) 

V 1/ 

Jq(z) = / + z-°e"(^2-^i)(^)E2i, -2 G Af nL>(0,r), (5.6) 

3. As z — > 0, Q has the same behavior as T has, see item 4. in the RH problem 3.2. 



25 



As we will see in Section 7, the radius ro will actually depend on n, namely ro = n~^/^, 
so that the parametrix Q will be defined in a disk shrinking neighborhood as n — t- oo. 

Note that we did not state a matching condition for Q. Usually one asks for a matching 
condition of the type 

Q{z) = Na{z){I + 0{l/n'')), 

as n — >■ oo, uniformly for z on the circle \z\ = vq, with some k > 0. In the present situation 
we arc not able to get such a matching condition. We will only be able to match Q{z) 
with Na{z) up to a bounded factor. Hence, it will be necessary to introduce an additional 
transformation, defined globally in the complex plane, as a last step of the Riemann-Hilbert 
analysis, see Section 7. 



5.3 Reduction to constant jumps 

The jump condition (5.2) can be reduced to a condition with constant jump matrices as 
follows. We put 

K{z) = diag(l,2", 1) diag (^e'^^i^^), e-"^^^^), e-"^^^^ ^ (5.7) 

with = |2|"e*"^''s(^) and argz G (0, 27r) is defined with the branch cut [0, +00). Then the 
jump matrices Jq from (5.3)-(5.6) factorize as 

jQ(z) = A-i(^)J^(^)An,+ (^), (5.8) 

where 

/I \ 

4(x) = -e-«- , X G A2, (5.9) 
\0 e-"'^' / 

j3(z) = 7 + e±"-E23, ^eAf, (5.10) 

/O 1 0\ 

J§(x)= -1 , xeAi, (5.11) 

\o 1/ 

J%{z)=I + E2i, zeAf, (5.12) 



Observe that the jump matrices Jq agree with the jump matrices J$ in the RH problem 1.3 
for see Figure 3, except that the jump matrices J$ are on six infinite rays S$ emanating 
from the origin. We will therefore look for Q in the form 

Qiz) = Eniz)^Mnizy, r„(z))A„(z) (5.13) 

where En{z) and r„(z) are analytic in D(0, ro) and where /„ is a conformal map on D{0,ro) 
that maps the contours n -D(0, ro) into the six rays S$ so that [0, ro) is mapped into the 
positive real axis. 

For any choice of conformal /„, and analytic En and r„, the matrix valued Q defined by 
(5.13) will then satisfy the required jumps (5.2). 
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Figure 12: The contours of integration Tj,j = l,... 4, used in the definition of the functions 
Pj. 

5.4 Definition of $0,(2;; r) 

We next construct the matrix valued function ^a{z) = $0(2:; r) that solves the RH problem 
1.3 stated in the introduction. 

As already mentioned in the introduction we use the following third order linear differential 
equation 

zp'" + ap" -Tp' -p = 0, (5.14) 

with a > — 1 and r G C. Then 2; = is a regular singular point of this ODE with Probenius 
indices 0, 1, and — a + 2. There are two linearly independent entire solutions, and one solution 
that branches at the origin. 

Due to the special form of the ODE (5.14) (the coefficients are at most linear in z), it can 
be solved with Laplace transforms. We find solutions with an integral representation 

p(2) = C^t"-W/(2*^)e^*di, 

where T is an appropriate contour in the complex t-plane and C is a constant. A basis of 
solutions of (5.14) can be chosen by selecting different contours T. We will make use of four 
contours Tj, j = 1,2,3,4, defined as follows, see also Figure 12: 

(a) We let Fi be a simple closed contour passing through the origin, but otherwise lying 
in the right half-plane and which is tangent to the imaginary axis. Fi is oriented 
counterclockwise and we put 

pi{z)= f r-W/(2*')e^*dt. (5.15) 

We choose the branch of t""^ = |^|a-3gi(Q-3)argt ^j^j^ _^/2 < argt < 7r/2. 

(b) F2 is the refiection of Fi in the imaginary axis, oriented clockwise. We put 

P2(z) = e-"'^^ /" t"-V/*e^/(2*')e^*dt. (5.16) 

In (5.16) we define the branch of with 7r/2 < avgt < 3ir/2. 



27 



(c) Fs is an unbounded contour in the upper half-plane that starts at infinity at an angle 
where Re (zt) < as f — >■ oo, and ends at the origin along the positive imaginary axis. 



In (5.17) we take with < argt < tt. 

The condition Re (zt) < is necessary to have convergence of the integral. This 
condition can be met with a contour that is in the upper half-plane if and only if 
— 7r/2 < argz < 37r/2. Therefore ^93 is well-defined and analytic in C\i]R_, that is, with 
a branch cut along the negative imaginary axis. 

(d) r4 is similar to T3 , but in the lower half-plane. It is an unbounded contour in the lower 
half- plane starting at infinity at an angle where Re {zt) < as t — >■ 00, and it ends at 
the origin along the negative imaginary axis. We put 



In (5.18) the branch of ^ is defined with — vr < argt < 0. 

Then is well-defined and analytic in C\iM+, thus with a branch cut along the positive 
imaginary axis. 

With these definitions it is clear that pi and p2 are entire functions, while ps and p4 have 
a branch point at the origin. The four solutions are not linearly independent, but any three 
of them are. 

We define in each of the sectors determined by the six rays: arg z = 0, ±7r/4, ib37r/4, tt, 
as shown in Figure 3, as a Wronskian matrix using three of the functions pj. 

Definition 5.1. We define $q. in the six sectors as follows. 



We put 




(5.17) 




(5.18) 






< arg 2; < 7r/4, 






(5.24) 



(5.23) 



(5.21) 



(5.22) 



(5.20) 



(5.19) 



28 



The scalar factor is needed in (5.19)-(5.24) in order to have the exact asymptotic 
behavior (1.22) in the RH problem 1.3. 

The functions pj clearly depend on a and r, although we did not emphasize it in the 
notation. 

Proposition 5.2. Let a > —1 and r G C. Then $0,(2;; r) as defined above satisfies the RH 
problem 1.3 stated in the introduction. 

Proof. It is a tedious but straightforward check based on the integral representations (5.15)- 
(5.18) that $Q has the constant jumps on six rays in the complex z-plane as given in (1.20). 

The asymptotic properties (1.22)-(1.23) follow from a steepest descent analysis for the 
integrals defining the functions pj. We give more details about this in the next subsection, 
where we also describe the next term in the asymptotic expansions of (1.22)-(1.23) since we 
will need this later on. 

The behavior (1.24)-(1.26) at follows from the behavior of the solutions pj of the ODE 
(5.14) at 0. Since p\ and p2 are entire solutions, they are bounded at 0. The solutions ^3 and 
P4 satisfy 

p,(z) = o(z2-"), p;.(z) = o(zi-"), p;'(z) = o(z-"), i = 3,4, 

as 2; ^ 0, which can be found by analyzing the integral representations (5.17) and (5.18). 
This proves (1.24)-(1.26) in view of the definition of in in Definition 5.1. □ 



5.5 Asymptotics of 

As before we define uj = e^'^*/^ and 9k as in (1.21). We also put 

('ga7ri/3 



La{z) 



-a/3 



,^1/3 

1 









OJ 








1 






u 







z-^/^j 





— OJ 


1\ /e-"'^*/^ 







-1 








-c^2 


; I 






(5.25) 



for Im 2; < 0, 

where all fractional powers are defined with a branch cut along the negative real axis. Define 
also the constant matrices 



M+(t) 



r(r2 + 9q; - 9) 
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diag {uP'^u, 1) 



+ 



IT 



:diag(|a;-"/^a;"/^l) 

X ■ V ^ + 9q;-9) . 2 -,\ 
Ma (r) = ^ diag {u, u:\ l) 



-uj 1 

w2 -1 I diag l) , 

^ ' ^ ^ 



(5.26) 
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'0 -a;2 -r 

+ ^diag (tj"/2^a;-"/^ \u -1 diag l) 



(5.27) 



kuj uj^ 
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Lemma 5.3. Let a > — 1 and r G C. Then we have, as z ^ oo, 

(gOiiz) \ 

e^^W (5.28) 
e^3(.)y 

for Im 2: > 0, 




























^^{z-t) = ^L^{z)\^I+^^ + 0{z-''')]\ e^^W I (5.29) 

for Im 2: < 0. 

The expansions (5.28) and (5.29) are valid uniformly for t in a compact subset of the complex 
plane. 

Proof. Wc apply the classical steepest descent analysis to the integral representations (5.15)- 
(5.18) of the functions pj. We set (j{t;z,T) = zt + rjt + l/(2i^). The saddle points are 
solutions of 

do T 1 , , 

^=.-^-p=0. (5.30) 

As z — )• 00, while r remains bounded, the three solutions to (5.30) have the following expan- 
sion: 

t^ = t^(^;r)=a;2fe^-V3 + ,^fcI^-2/3 + 0(^-4/3)^ A; = 1,2, 3, (5.31) 
and the corresponding values at the saddles are 

a{tk{z; t)-z, t) = \u?''z^l^ + rJ'z^l^ - ^ + 0{z-^l^) 

2 

= Qk{z;T)-'^— + 0{z-^l'^), asz^oo, (5.32) 

with introduced in (1.21). 

If Cfc is the steepest descent path through the saddle point i^, we obtain from (5.32) and 
standard steepest descent arguments that 



/ t"-3e-/*eV(2*^)e^*dt = ± . / t"-3g--V6ge,(.;r)(i + o(z-V3)) (5.33) 

JCk ]l -%i{tk;z,T) 

as z ^ 00, where the it sign depends on the orientation of the steepest descent path. Plugging 
(5.31) into (5.33), and using the fact that 

3 f_, 2 \ 
^ + i^rt], 



dt^ 3 

we obtain as z ^ 00 



J t"-W/(2*')e^*dt = ±y^e--'/6(u;*^2-V3)" ' e«fc(-;-) (l + 0(z-V3)) . (5.34) 
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Figure 13: Deformation of the contours of integration Fi, F3 and F4 into the steepest descent 
paths for the integrals defining pi, and p^ in the case where < argz < 7r/4. 



Consider now ^aiz) as defined in (5.19) for < argz < tt/4. We see that only pi, p^ and 
Pi play a role in this sector. The corresponding contours of integration can be deformed into 
the steepest descent paths through one of the saddle points as shown in Figurel3. 

Hence, in the given sector, the functions p\ , p^ and p4 have an asymptotic behavior of the 
form (5.34) for some particular k and some choice of the it sign, and multiplied by e^"'^* in 
case of p3 and pi, see the formulas (5.15), (5.17), and (5.18). This will lead to the asymptotic 
expansion for the first row of (5.28) for < z < 7r/4, except for the determination of M+(r). 
The second and rows can be dealt with similarly. Here we have to consider the first and 
second derivatives of pi, ^3, and p4, which have similar integral representations (5.15), (5.17), 
and (5.18), but with a replaced by a + 1 and a + 2, respectively. 

The other sectors can be analyzed in a similar way. Tracing the behavior of the dominant 
saddle point we find that the asymptotic expression just obtained remains valid in the full 
upper half plane, while in the lower half plane we find (5.29), again up to the determination 
ofM-(r). 

What remains is to obtain the constants M^(r) that appear in the 0{z-^/^) term in 
the expansions (5.28) and (5.29). This can be done by calculating the next terms in the 
asymptotic expansion of the integrals. Alternatively, we can use the fact that solves the 
first-order matrix-valued ODE 



Substituting into this the asymptotic expansions for $a and equating terms on both sides, 
we find after lengthy calculations (that were actually performed with the help of Maple) the 
formulas for M^(r) as given in (5.26) and (5.27). □ 

5.6 Definition and properties of f{z) and t{z) 

We will take the local parametrix Q in the form, see also (5.13) 

Q{z) = En{z)^a{n^^'^f{z)-'n}/'^T{z))\ z'^(,-n\2{z) esfe, (5.35) 

e-n\3{z)) 
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where En is an analytic prefactor, f(z) is a conformal map defined in a neighborhood of 2: = 
and t{z) is analytic in z. Assuming that / maps the contour T,t H D(0, r) to the six rays S$ 
such that f{z) is positive for positive real z, then it follows from the above construction that 
Q will satisfy the required jump condition. 

We are going to take f{z) and t{z) in such a way that the exponential factors in (5.35) 
are cancelled. That is, we want analytic f{z) and t{z) such that 

ek{f{z)-T{z)) = \k{z) - ^^^—^y k = l,2, 3, (5.36) 
for Imz > 0, while for Imz < 0, 

eiif{z);T{z)) = \2{z) - ^^^^ - 27rz (5.37) 

^2(/(^);r(z)) = Ai(z) - ^^^^ + 27ri (5.38) 

Hf{^y,r{z)) = Xsiz) - ^^^—y (5.39) 

where the functions 9k were defined in (1.21). 

To define f{z) and t{z) we use the functions h{z\t) and g^{z\t) from Lemma 2.8. 

Definition 5.4. We put 

f{z) = f{z; t) = z[Mz; t)f', r{z) = t{z; t) = (5.40) 

where as usual we take the principal branches of the fractional powers. We write f{z] t) and 
g{z;t) in order to emphasize their dependence on the parameter t. 

Lemma 5.5. There exist vq > and S > such that for each t E {t* — S,t* + S) we have 

that z ^ fiz;t) is a conformal mapping on the disk D{0,ro) and z 1— t- T{z;t) is analytic on 
D{0, ro). The map z 1— )• f(z; t) is positive for positive real z and negative for negative real z. 
In addition, we have 

T{z;t) = 0{t-t*) + 0{z) ast^t* andz^O. (5.41) 

Proof. Because of (2.26) we have that /3(0) > and so f{z) defined by (5.40) is indeed a 
conformal map in a neigborhood of z = which is positive for positive values of z. Also t(z) 
is analytic in a neighborhood of 2; = 0. 
We have 

f3{z;t) = f,{z,t*) + 0{t-t*), g3{z;t)=gs{z,t*) + 0{t-t*), ast^f, (5.42) 

uniformly for z in a neighborhood of 0, and 

f3iz;t*) = {c*)-''/^ + 0{z), gsiz-X) = 0{z) as z ^ 0. (5.43) 

This follows from the definitions (2.14) of c and p, equation (2.16), and the definitions of fj 
and gj, j = 1,2,3. Expansions (5.43) also use (2.26) and the fact that p = when t = t*. 
Then (5.41) is a consequence of the previous expansions and the definitions of f{z) and t{z) 
in (5.40). □ 
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5.7 Definition and properties of the prefactor En{z) 

The prefactor En{z) in the definition of Q{z) in (5.35) should be analytic and chosen so that 
Q is close to A^^^ on \z\ = nT^/"^. In view of the expansion of $ given in Lemma 5.3, wc set the 
following definition. We use ro > as given by Lemma 5.5 and we assume i G (i* — 5, t* +8). 

Definition 5.6. We define for z G r>(0,ro) \M, 



/I 0\ 

En{z) = -iVSNaiz) z-^ L-\n^/^f{z)), 
\0 ij 

where La has been introduced in (5.25), and Na is described in Section 4. 
Lemma 5.7. and E~^ have an analytic continuation to D{(},ro). 
Proof. Taking into account (4.1)-(4.2) we see that for Nq, defined in (4.13), 



(5.44) 



1 0> 

Na,+ {x) = Na,-{x) 1-10 
1, 



for X G (0,ro) 



and 



10 

Na,+ix) = Na,-ix) \ -e-"-^ 



for X G (— ro, 0). 

For La we find the same jump matrices. Indeed, for x > 0, we have by the definition of 

La-{x)''^La,+ {x) = 




and for x < 0, 
La-ix)~^La,+ {x) 











e 






^ 






0\ 




H 





1 






{0 








^—aivi/S Q 








1 -1 




LO —Li? 












V 





33 



jaTTi/s 0\ /I \ /e"''*/^ 0^ 

^g-2a7ri/3| q g-a^S Q -1 e'^^^/^ q 




l/V0 10/\0 1 



10 
= 10 -e-"'^* 
^0 6-°^'^^ 

Thus, the jumps for Na and La are the same. Since / is a conformal map on D(0,ro) with 
f{x) > for X G (0,ro) and /(x) < for x < 0, it follows that Na{z)L-^{n^/'^f{z)) has 
an analytic continuation to D(0, ro) \ {0}. As a result, we conclude that has an analytic 
continuation to D(0, ro) \ {0}. 

We show that the isolated singularity at the origin is removable. Indeed, by Lemma 4.3 
we have 

Na{z) = as z ^ 0, 

and by the definition of 

z-'^/^L-\z) = 0{z-^/^) asz^O. 

Thus 

a/3 



Z 



--''Na{z)L-\n'l'f{z)) = Na{z) {n^/^ f{z)-'^/'L-\n'/' f{z))) 

= 0{z-^/^) asz^O. 

It follows that the singularity of the left hand side at 2; = is removable and thus En{z) is 
analytic in D{0,ro). 

Recall that by Lemma 4.3 we have that det Na = dot Na = 1- From (5.25) we get that 



det La{z) = SiVSz 

Thus by (5.44) 

det En{z) 



n' 



z 



which is analytic and non-zero in a neighborhood of z = 0. Thus E~^[z) is analytic in the 
neighborhood D{0,ro) as well. 

This completes the proof of the lemma. □ 

Having defined f{z), t{z) and En{z) we then define the local parametrix Q as in formula 
(5.35). 



6 Fourth transformation of the RH problem 

In the next transformation we define 



S{z) 



T{z)P{z)-\ for zeD{q*,rg), 
T{z)Q{z)-^, for z G D(0, 71'^^), 
T{z)N~^{z), elsewhere. 



(6.1) 
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Figure 14: Jump contour S5 for the RH problem for S. The disk around is shrinking with 
radius n~^/^ as n ^ 00. The contour S5 is also the contour for the RH problem for R. 



where we use the matrix- valued functions Na from (4.9), P constructed in Section 5 in the 
fixed neighborhood D{q*,r) of q*, and Q given by (5.35) in the shrinking neighborhood 
D{0, n~^/^) of the origin. 

By construction, S{z) is piece- wise analytic and has jumps across the contour Eg shown in 
Figure 14, with a possible isolated singularities at and q*. The singularity at q* is removable 
which follows from the properties of the Airy parametrix. We now check that the origin is 
not a singularity of S{z). 

Lemma 6.1. The singularity of S at z = is removable. 

Proof. We give the proof for the case a > 0. 

Consider z ^ with Imz > and outside the lenses around A2 and Ai. By the RH 
problem 3.2 for T, we have that T remains bounded there. We show that remains 
bounded as well. 

By (5.35), we have 



By Lemma 5.7 we know that En{z) ^ is analytic and thus bounded as z — >■ 0. Also the 
functions A,- are bounded as z ^ 0. Also 



is bounded as z ^ in the region under consideration because of the condition 4. in the RH 
problem for sec (1.25) and the fact that det = z~°', see (1.27). 

Wee conclude that remains bounded as 2; —> in the region in the upper half-plane 
outside of the lenses. 






The other regions can be treated in a similar way and the lemma follows. 



□ 



We find the following RH problem for S 
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RH problem 6.2. 1. S is analytic outside the contour S5 shown in Figure 14. 



2. On S5 there is a jump relation 

S+iz) = S.iz)Jsiz) (6.2) 

with jump matrix Js{z) given by 

Js{z) = N^{z)p-\z), for \z-q*\= r^, (6.3) 

J5(^) = Na{z)Q-\z), for |z| = ro = n-^/^ (6.4) 

Js{z) = Na{z)JT{z)N~^{z), elsewhere on E5 (6.5) 



3. S{z) = 1 + C(z~^) as z ^- 00. 

Recall now that we are interested in the limit (2.1) where n ^ 00, t ^ t* and c*r = 
v}/'^{t* — t) remains fixed. The jump matrices J5 in (6.3)-(6.5) depend on n and t, and we 
would like that they tend to the identity matrix in the double scaling limit (2.1). This turns 
out to be the case for the jump matrices (6.3) and (6.5). However, this is not the case for 
(6.4) as will be shown later on. 

We start with the good jumps. 

Lemma 6.3. In the limit (2.1) we have 

Js{z) = 1 + {n~^) uniformly for \z — q* \ = rq, (6-6) 
and for some c> ( depending on a ), 

Js{z) = I + O \ \ -\-^ I j uniformly for z G S5 outside of the two circles. (6-7) 

Proof. The behavior (6.6) follows from (6.3) and the matching condition (5.1) for P. 

In view of Lemma 2.9 and the asymptotic behavior of the A functions we find that for 
some c > 0, 

Re(A3- A2)(z) > c|z|^/2^ \i:»(0,l), (6.8) 

Re(Ai-A2)(z)>c, A±\(L>(0,l)UD(g*,rq)) (6.9) 

Re(A2 - Ai)(z) > c|z|, ze{q* + rq,oo). (6.10) 

According to (3.5), (6.5) and (6.8), for z G \ D(0, 1), 

Js{z) -1 = e±"'^^z-"e'^(^2-A3)(^)jv^(^-)^23iV-i(z) = O (l^r^e"'^"!^!'^') 
for some c > 0. 

Analogous considerations on the lips A]*" \ (-D(0, 1) U D{q*,rq)) and on {q* + rq, 00), ap- 
pealing to formulas (3.7)-(3.8) and (6.9)-(6.10), show that there exists some c > 0, such 
that 

Js{z) =1 + (|z|He— 1^1''') , ze^s\ iD{0, 1) U D{q*,rq)). 
What remains is to estimate Jsiz) on the lips of the lenses near for n~^/^ < \z\ < 1. 
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For t = t*,we obtain from Lemma 2.8 that there exists a constant ci > such that (recall 
that p = when t = t*) 

Re (A3 - \2){z; f) > Cl|zp/^ zeAfn D{0, 1), 

Re (Ai - X2){z; f) > ci\zf/^, z e Af n D{0, 1), 

Moreover, (2.25) and (5.42) imply that Xj{z,t) = Xj{z;t*) + z^/^ 0{t - t*) as t t*. Thus, 

Re(A3 - A2)(z;t) > ci\z\'^/^ - C2\z\'^^^\t - t*\, z G A^nD(0,l), 

Re(Ai - A2)(2;t) > cilzl"^/^ - C2\z\^/^\t - t*\, z G Af nL>(0,l), 

for some C2 > 0. Since t — t* = 0(n~^/^) we conclude that 

Re (A3 - X2){z; t) > c^n-^/^, z e Af D D{0, 1), n^^/^ < \z\ < 1, 

Re (Ai - X2){z; t) > cm"^/^, z e A"^ n D{0, 1), n~^/^ < \z\ < 1, 

for some positive constant C3 > and n large enough. 

Now, for zeAfn D{0, 1), \z\ > n~^/^, and using (4.13), we get 

Jsiz) -1 = e±"-^e'^(^2-A3)(.)^^(^)^23iV-i(^). 
By Lemma 4.3, Na{z) = O {\z\-^/^), N-^{z) = O {\z\-^/^) as \z\ 0, so that 

J5(z)-/ = 0(e— 

for some c > 0. Analogous conclusion is obtained on n D{0, 1), n~^/^ < \z\ < 1. 

Gathering both estimates and replacing them by a weaker uniform bound we obtain 
(6.7). □ 

We next analyze the jump matrix (6.4) on \z\ = n~^/^ again in the double scaling limit 
(2.1) 

We have = NaQ ^ where Na is given by (4.9) and Q is given by (5.35). All notions 
that appear in these formulas depend on i or n (or both). For example, A^^ depends on t 
since the endpoint q is varying with t, and tends to q* as t ^ t*. Indeed, q = q* + 0{t — t*). 
Also the matrix Cq from (4.9) and the constants Ki, K2, K3 from (4.7) depend on t and tend 
to limiting values corresponding to the value t* at the same rate of 0{t — t*) = 0{n~^/'^). We 
denote the limiting values with *: 

Ca = C*+0(t-r), Kj = K*+0{t-t*), j = 1,2,3, (6.11) 

and these quantities appear in the formula (6.12) below. 

Proposition 6.4. In the limit (2.1) we have that 

j^(^) = / _ ^h^Ml + 0(n-i/6) (6.12) 

uniformly for \z\ = n~^/^ , where 

/K*A /I 

Ml = C*jK*,\{Kt Kl i^l) -i\{C*^r' (6.13) 
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and 

hn{z;t)=T{z;t) ^ '^^^ . (6.14) 

Proof. The local parametrix Q from (5.35) depends on both t and n. The functions Xj come 
from the Riemann surface and therefore depend on t, but only in a mild way. Of more 
importance is the dependence of the functions f{z) = f{z;t) and t{z) = T{z;t) on t, see 
Lemma 5.5. 

From (5.41) we have that 

nVV(z;t) = 0(1), (6.15) 

in the double scaling limit (2.1) uniformly for \z\ = n~^/^. [This is in fact the reason why we 
need the shrinking disk of radius n~^/^.] 

We also have that n^^'^ f{z;t) remains bounded as \z\ = n~^^'^, However, rfi/"^ f{z;t) is 
growing in absolute value and is of order n uniformly for \z\= n~^/^. Therefore we can apply 
the asymptotic formulas from Lemma 5.3 to $q, {n^/'^f{z); n^/'^T^z)) and we find for Imz > 0, 



(n^/V(^);nVV(z)) 

X ^e,in^/'f{zyy/^r{z)) I (6.16) 

^ ge3(n3/2/W;ni/V(2))^ 

uniformly for \z\ = n~^/^, where we recall that 9k{z;T), k = 1,2,3, also depend on r, see 
(1.21). By (1.21) and (5.36) we actually have 

9, {n^/'f{zy,n'/'T{z)) = ne,(/(z); t(z)) = nXk{z) - 

for A; = 1, 2, 3, and Imz > 0, by our choice of f{z) and t{z). 
Thus by (5.35) and (2.1) we have 

Q(.) = E.(.)-j=L„(n3/V(.))(/+^^g|^^ (I r 

and then by inserting the definition (5.44) of En{z), we obtain 

QU) = iV,(.) (I A oj (/ + + 0(„-/3)) g ^ oj . (6.17) 

uniformly for \z\ = n~^/^ with Imz > 0. For \z\ = n'^l'^ with Imz < we obtain the same 
formula (6.17) but with M+ replaced by M~. 
Then for \z\ = ra~^/^, 

Js\z) = Q{z)N-\z) = I + N^[z)^^j:^j^^N^\z) + 0{n-''^) (6.18) 
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where we recall the definition (4.13) of N^- The entries of Na{z) and its inverse are of order 
= n^/^ by (4.14). Therefore the error term has gone up from 0(n~^/^) in (6.17) to 
0(n-i/3) in (6.18). 

In the evaluation of Naiz)M^{n^/'^T{z))N^^{z) we encounter the following matrix 

Ma = Ca\K2 \ {Ki K2 K^) [ -z (6.19) 
\KsJ \0 -i 0/ 

which is a 3 X 3 rank one matrix depending on t, but not on z. The matrix Aia is in fact 
nilpotent, 

Ml = 0, (6.20) 

which follows from (6.19) and the property (4.16) of the numbers Kj. Prom (6.11) and (6.13) 
we also find that 

Ma = Ml + 0{t~t*) ast^t*. (6.21) 
We first prove the following lemma. 
Lemma 6.5. We have for Imz > 








0' 



















Nc,{z)\ 00 N-\z) = Sc-'^'-'MaZ-'^'' + O z"^/^ (6.22) 



and 






—UJ 














-:■) 


1 : 




-a;2 












0^ 

Na{z) I 

= -3V3zc-^/=^A^«z-2/3 + O (6.23) 

OS z — 0. 

Proo/. We obtain (6.22) from (4.9), (4.12), Lemma 4.4, and the fact that 



3. 



Prom (4.9), we obtain that the left-hand side of (6.23) is equal to 

fif{z) W ^ -a; 1 \ /er"(^) 
V i^{z)l CO Oj \ ^^"(z), 

with functions 

^i(^)=a;-i/2gGi(.)^ 6(^)=^'/'^-'e^^(^\ e3(-^) = e^^^^^, 















CO 
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with Gj{z) = rj{wj{z)) defined in (4.10)-(4.11). Using these expressions and Lemma 2.5, we 
find the remarkable fact that 



Then, using (4.12) and the fact that 



1 + O (z^'^^ , z^Q. 






—OJ 









-:■) 


-c^2 







[u"^ u l)\ uj"^ -1 a; = -SVSi, 
a; y \ 1 / 

we obtain (6.23) □ 

Wc continue with the proof of Proposition 6.4. 

From tlie lemma and the formula (5.26) for , we obtain 

N^{z)M+{n"\{z))N-\z) 

= ,i/2 r(.)(nr^(z) + ^ ^ (^-i/3-| (6_24) 

Inserting this into (6.18) we get for \z\ = n~^l'^ and Imz > 

A similar analysis for Imz < will show that the same formula (6.25) also holds for Imz < 0. 

Note that for \z\ = rT^I"^ we have ^^/i^^^^^i/a = 0{n~^l^\ where we use that j{z) is a 
conformal map. In fact, by (2.26) and (5.40) 

j{z- 1) = /'(O; t)z + 0{z') = {(c*)-' + Oit -n)z + Oiz^) 

which implies that in the double scaling limit (2.1) 

fiz;t) = ic*)-'z + Oin-') 

for \z\ = n~^/^. Using also (6.20) we then obtain from (6.25) that 

M.) ^ I - ^'^";;'^g3^ '"^ M^.-^'^ + 0(n-/°), 1.1 = (6.26) 
which implies (6.12) in view of (6.21). This completes the proof of the Proposition 6.4. □ 

Now recall that in the double scaling limit (2.1) we have that n^/^r(z; t) remains bounded 
for \z\ = n"^/^, see also (6.15). Then also 

T{z;t)inT\z;t) + 9a) ^ ^ 

c*z ^ ' 

and it follows that the term ^^^^^X* in (6.12) remains bounded for |2;| = n but does 
not tend to as n — oo. Therefore the jump matrix J5 on \z\ = n~^/^ does not tend to the 
identity matrix as n — >^ 00. 
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7 Final transformation 



We need one more transformation. What will help us in the final transformation is the 
identity (6.20) for M.a, which also holds for the limiting value, namely {M.*^)^ = 0. The final 
transformation S R is similar to the one in [16] and is defined as follows. 

Definition 7.1. With the notation (6.14) we define for z G C \ S5, 

R{z) = S{z) ( I - M*^ , \z\ > n-V2, (7.I) 



Note that the transformation S R is & global transformation which modifies S in 
every part of the complex plane. Then R satisfies the following RH problem on the contour 
Efl = E5, see Figure 14. 

RH problem 7.2. 1. R is defined and analytic in C \ YIr. 
2. On we have the jump 

R+ = R^Jr (7.3) 

with 

Mz) = {i+ ^ m:) js{z) [i + hn{z;t)-Km ^ ^7.4) 

for 1^1 = n~^/^, and 

Mz)=(l+'^M-yJsi.z)il-'^M'y, (7.5) 



elsewhere on T^r. 
3. R{z) = 1 + 0{l/z) as z ^ 00. 

All properties in the RH problem 7.2 follow easily fom the RH problem 6.2 for S and the 
definition (7.1)-(7.2). For (7.4) and (7.5) one also uses (6.20) and (6.21) which imply that 
for every constant 7, 

{i-^Mir' = i+^Mi. 

Under the transformation S ^ R the jumps on the part of T^r outside of the circle 
\z\ = n~^l'^ are not essentially affected. We have the same estimates as in Lemma 6.3: 

Lemma 7.3. In the limit (2.1) wt have 

Jr{z) = 1 + (n~^) uniformly for \z — q*\ = Vq. (7-6) 

and 

Jr{z) = I + O ( ^ ^ 1 uniformly for z G T^r outside of the two circles. (7.7) 



41 



Proof. By (6.6) and (7.5), for l-z - g*| = rg, 



Jr={I + 



/+ / + 



Ml^ {l + 0{n-'))(^I-^^Ml^ 
= I + 0{n-^) . 

Analogous calculations yield (7.7). In both cases the fact that (A4* )^ = is crucial. □ 
Lemma 7.4. In the limit (2.1) we have 

Jji{z) = I + (n"^/^) uniformly for \z\ = n'^/^. (7.8) 

Proof. Using (7.4) and (6.26) and the fact that (A^^)^ = 0, direct calculation yields 

Mz) = /+(/+ ^ Ml^ O (n-V6) {l + K{z;t)-Km j^^^ 



and (7.8) follows. □ 

As a result of the estimates on Jr we may now conclude that in the double scaling limit 
(2.1) 

^W=-^ + o( „Vaa + |,|) ) P-S) 

uniformly for ^ G C \ Yjr. See [2, Appendix A] for arguments that justify this, also in a 
situation of varying contours. 

8 The limiting kernel 

8.1 Expression for the critical kernel 

We start from (1.12), which gives the correlation kernel Kn{x,y;t) in terms of the solution 
of the RH problem for Y. Following the transformation Y t-^ X t-^ U i-^ T, we find that for 
x,y > and x,y E {0,q), 



Knix,y;t) = ^^.^^ _ (-e-'*^i.+(2') y«e-"^2.+(2/) g) T-\y)T+{x) | a;-"e"^2,+(a;) 



oraAi_+(x-) 





(8.1) 

For z inside the disk of radius n we have by (6.1), (5.35), (5.7), and (7.2), that 

T{z) = Siz)Q{z) 

^ fj _ T(z;0(nT^(z;t) + 9«)-r(0;t)(nT^(0;t) + 9a) A 



X 



En{z)^a{n^^^f{z; t); n'/\{z; i))A„(z)e3*("-*) , 
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Thus, i{ < x,y < n we get by plugging (8.2) into (8.1), 



Kn{x,y;t) 



2n(x — y) 
g 3t(l-t) 

27ri(x — y) 



1 1 o)<i>-V(n3/V(y;t);n^/My;i)) 



X Ri.) (l - ^(^; t) + 9a) - r(0; t) (nr^(0; t) + 9a) ^j ^^^^^ 



X $c«,+(n2/V(^; 0; "'/'^(a^; *)) 1 1 1 (8-3) 



which is an exact formula. 

Now we take the double scaling limit n — > oo, t — )■ t* such that 

c*T = n}^^{t* — t) remains fixed. 

We also replace x and y in (8.3) by 



c X 



Xn. — 



n' 



3/2' 



- „3/2 



(8.4) 



(8.5) 



with a;, y > fixed. For n large enough we then have that Xn and y„ are less than n ^1'^ . We 
study how the various factors in (8.3) behave in this limit. 

Lemma 8.1. Let x,y > he fixed. Then we have in the double scaling limit (8.4) with Xn and 
yn given by (8.5) 



n 



3/2/(xn; t) = x{l + 0(n-i/2)), n3/2/(yn; t) = y{l + 0(n-V2)), 



(8.6) 



and 



n 



1/2 



T{xn; t)=T + 0{n-'/^), n'/^yn, t)=T + 0{n-'/^) 



(8.7) 



as n oo . 

Proof. By (5.40), (5.42) and (5.43) we have 

f{z; t) = z[Mz; t)f^ = z [{c*)-'^' + 0{z) + 0{t - t*)] ^'^ = ^[l + 0{z) + 0{t - f)] 

as z — 7- and t t* . This readily implies (8.6). 
Again by (5.40), (5.42) and (5.43) we have 



Tiz-t) 



93{z;t) 

hiz-tyi^ 



e f^'g^iz-t) [l + 0{z) + 0{t-f)] 



as z — and t ^ t* . Then from the definitions in Lemmas 2.5, 2.6, and 2.8 it is not difHcult 
to verify that g3{z;t) is analytic in both arguments with 53(0, t*) = 0, so that 



93iz,t) 



dgs 

at 



(o,r 



*\2 



{t-t*) + Oiz) + 0{t-t*) 
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as 2 — >■ and t t*. By (2.26) we have 53(0; t) = 3pc and using the dependence of p and 
c on t as given in (2.14) we find, after some calculations, that 

^(o,f) = -(c*)-^/^ 

Hence r(z;t) = (c*)~^/^ {t* -t) + 0{z) + 0{{t - t*)^) and 

ni/V(x„; t) = (c*)~^ (f _ t) + 0{n-^) + n^'^0{{t - t*f) 
which by (8.4) indeed leads to (8.7). □ 
Lemma 8.2. Under the same assumptions as in Lemma 8.1 

\ 9c* y 

X -R(X„) I / ^-^^j Maj En{Xn) I- (8.8) 

Proof. For 2; = 0{n~^/^) we have by Cauchy's theorem and (7.9) 

R'(z) = ^ I ^^'^^ ~,l ds = C»(nV3) as n ^ 00. 

27rz J|^l=„-i/2 [s-zY 

Therefore i?(x„) - R{yn) = 0((x„ - 2/n)«^/^) = 0(n-^/^) and so 

i2-^(y„)i?(x„) = / + i?-^(y„)(i?(x„)-i?(y„)) = / + C'(n-^/6) as n ^ 00. (8.9) 

where we use that R~^{y n) remains bounded as n — >■ 00, which also follows from (7.9). 
Let us write 

_ Tjxn; t)inT'^{xn; t) + 9a) - r(0; t)(nT^(0; t) + 9a) 

Pn{Xn,t) .— {Q.W) 

and similarly for Pn{yn,t). Then explicit calculations (done with the help of Maple) show 
that 

+ 3a T^(a — 2) 1/9 1. /^,,s 

p4x„, t) = - s4(„Vl)(c-)^ " ' + ' 

and similarly for Pn{yn,t). Thus 

/9n(yn, i) - Pnixn, t) = 0{n~^) as n ^ DO. (8.12) 

Using (8.9), (8.11), (8.12) and the fact that {Mlf = 0> we see that (8.8) will follow from 
the following three estimates 

E-\yn)En{xn) = 1 + OirT^I'') (8.13) 

E-\y^)0(n-'l'^)Er,{xn) = 0(n-i/6) (8.14) 

E-\yn)MlEn{xn) = 0(v}l'') (8.15) 

as n — >■ DO. 
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By (5.44), the analytic factor En{z) depends on n mainly via the argument n^^"^ f{z] t) of 
L~^. By (4.13), (5.25), and (5.44), we can factor out the dependence on np/'^f{z;t), and we 
obtain 



where 



,g-Q7ri/3 0\ 1 w\ /z-1/3 

e"''*/^ w 1 1 

l/\lll/VO 



Fa{z) = T^Naiz) ( e'^'^^/^ 1 I w 1 a;2 I I 1 | (8.17) 



is analytic around z = and depends on n in a very mild way only, namely via the dependence 
of Na{z) on the endpoint q which is only slightly moving with n. 

The scalar factor in (8.16) will appear in the products (8.13)-(8.15) in the form 



q/3 



= l + 0{n-^l^) (8.18) 



where the estimate follows from (8.5) and (8.6). Thus by (8.16) 

1/3 ^ /„3/2r. -1/3 



^nHyn)£^n(x„) = (l + 0(n-V2))diag^!^!^!i^ 1 (l^^^^) 

xF-H2/n)Fa(Xn)diag((!l!^!|^)-'/' 1 [^^^^^f' 

3 2 1 /3 

where F-^{yn)Fa{xn) = l + 0{xn-yn) = l + 0{n-^/'^). Since the two entries ( "^^'^^^"'*^ ) 

and ^ " ^ f^^"'*^ j ^ in the diagonal matrices grow like C>(n^/^) we find (8.13), where we also 
use (8.18). 

We similarly have 



(nyy{y„;t) 













xF-i(y„)0(n-V6)F„(^„)diag(f!l!^!fe^)r'/' 1 (n^^lfi^V^' 

Since F~^(yn) and Fo,{xn) remain bounded as n ^ oo, and the two diagonal matrices are 
0{n}/^) we obtain the estimate (8.14). 

To prove the final estimate (8.15) we note that 



^n^(yn)>t:^n(x„) = (l + 0(n-V2))diag(^(l^!^!^ 1 (^^^^^)"'^' 

xF-\yn)MlFUxn)dmg((^!f!l^y'^' 1 (^V^!!^^^'^ (8.19) 
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which would lead to 0{n) as n ^ oo if we use the same estimates as above. However, by the 
form of the right hand side of (8.19), we see that it is only the (1, 3) entry E^^ {yn)M'^En{xn) 
that could grow like 0{n). The other entries are Oin^l'^^ as claimed in (8.15). 
We have by (8.19) and (8.17) 



' \ Xn / \ / ' 



O {\\N-\yn)M*Mxn)\\) 



where in the last step we used (8.6). Both matrices N~^{yn) and Na{xn) grow like n^/^, see 
(4.5). However 

M*o,Na{xn) = 0{1) as n ^ oo 
which follows from (4.13) and the fact that by (4.12), (6.19), 

/Ki\ /I 0\ /KA 

M^CaNoiz) = Ca\K2\ (i^i ^2 i^s) -i c'^/^ K2 {uj^ UJ l) z'^'^ + 0(1) 

\K^I \Q -i 0/ V^3/ 
= 0{l) as z ^ 

because of the relation (4.16) satisfied by the constants Kj. Thus N~^{yn)MaNa{xn) = 
0{n^/'^) and (8.15) follows. 

This completes the proof of Lemma 8.2. □ 

Prom (8.3), (8.5) and Lemmas 8.1 and 8.2, we obtain 

1* ^ T/^ i ^ ^ y J.* C T \ fy^j^j- , . 

where i^Q ** is given by (1.28). This proves Theorem 1.4. 
8.2 Proof of Theorem 1.2 

Let us analyze the expression (1.28) for the critical kernel. We define 

(P2 P3 Pi^ 
= P2 P3 Pi 
\p'i P'^ P'(j 

as an analytic matrix- valued function in C \ iM_. It is the analytic continuation of the 
restriction of to the upper sector 7r/4 < argz < 37r/4 to the cut plane C \ zM_. 
Then by the jump relations of $q,, see (1.20), we can rewrite (1.28) as 

KZ'\x,y;r) = ^-^^^{0 1 O) (^^i^')"^ (y)$TO | ] , (8.20) 
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for > 0. Clearly, 

^7{^) I I = I I . (8.21) 




The inverse of is built out of solutions of the differential equation 

xq" + (3 - a)q" - rq' + q = 0, (8.22) 

which is, up to a sign, the adjoint of the equation (5.14). 
Define the pairing 

\p{x)^(liy)] = {yQ"{y) - (« - 2)g'(y) - Tq{y))p{x) 

+ {-yq{y) + (« - l)q{y))p'ix) + yq{y)p"{x), (8.23) 

and denote [p, g](x) = \p{x),q{x)] which is the bilinear concomitant. Then 

4^b' oli^) — Pi^) {xq"'{x) + (3 — a)q"{x) — Tq'{x) + q{x)) 
ax 

+ q{x) {xp"'{x) + ap"{x) — Tp'{x) — p{x)) , 

which shows that if p and q satisfy the respective differential equations, then the bilinear 
concomitant [p, q] (x) is constant. 

To find the inverse of we need solutions, that we call qi, q2, and 53, dual to pi, p2, 
P3, satisfying 

[Pj,qk] = Sj,k, j,k = 1,2,3. (8.24) 
The inverse matrix is then given by 

I zq'^iz) - {a - 2)q'2{z) - Tq2{z) - zq'^iz) + {a - l)q2{z) zq2{zf 
($7)-^ {z) = zqUz) - (a - 2)q',{z) - Tq,{z) -zq',{z) + {a - l)q,{z) zq,{z) 
\zq'{{z) - {a - 2)q'^{z) - Tqi{z) - zq[{z) + {a - l)qi{z) zqi{z)j 

and 

(0 1 0) ($;if )-^ (y) = {y4{y) - (a - 2)g^(2/) - Tq^{y), -yq'^ + (a - l)g3(y), y©(y)) • 

(8.25) 

Hence, the solution ^3 is the relevant one for the critical kernel: by (8.20)-(8.25), we get 

-crit/ ,.^ _ \p2{x),q3{y)] 



K'"\x,y) = 



2m{x — y) 



Before continuing, let us build the dual functions for pj. The solutions qj{z), j = 1,2,3, 
of (8.22) admit integral representations 

qj(z) = Cj [ t-'*e-^/*e-^/(2*')e-^*dt. (8.26) 

When the variable z is positive, the contours Sj, j = 1,2,3, can be chosen as in Figure 15; 
note that the integrals converge since the contours T,j approach the origin tangcntially to the 
real axis and go to infinity along the positive real axis. We choose the main branch of t" in 
(8.26) with the cut in the t-plane along M+, and allow Si to go along the upper side of the 
cut. Observe that we may take the same branch cut in the definition of p2 in (5.16). 

In order to build an integral expression for the bilinear concomitant we need the following 
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Figure 15: The contours of integration Ej, j = 1,2,3, (in bold) used in the definition (8.26) 
of the functions qj{z), when z is positive. Also shown the contours Tj, j = 1,2,3 (dashed 
lines) . 

Lemma 8.3. Let p and q be solutions of (5.14) and (8.22) respectively, with integral expres- 
sions 

where T and S are one of the contours depicted in Figure 15, respectively, and 



V{t) = {a-2)\og{t) + -^+^. 



(8.27) 



Then 
(a) 



(b) 



j r^e^We^^dt = zp"{z) + (a - 1) p'{z) - Tp{z), 
j t-^e^^^h'^dt = -rzp"{z) + {z- r(a - 1)) p'{z) + (a - 2 + t'^)p{z), 
J t-^e-^^'^e-''dt = zq"{z) + (2 - a) q'{z) - Tq{z). 

\pix),q{y)] = 1 1 - ^^^^1^) e^W-^(^)e^*--^ dtds. (8.28) 

In particular, 

7r 7e t — s 

(c) If the only point of intersection ofV andT, is the origin, then \p,q]{x) = 0. 
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(d) // r and S intersect transversally at zq ^ 0, and if the contours are oriented so that F 
meets S in zq on the "— "-side of E, then 



Proof. Let us denote 



Hz) 



\p,q\{x) = 2m. 



Then h'{z) = p{z) = zp"'{z) + (a + 2)p"{z) — Tp'{z), where we have used the differential 
equation (5.14). Hence, h{z) = zp"{z) + {a + l)p'{z) — Tp{z) + c, and in order to find the 
constant c we compute 

h{z)-{zp"{z) + (a + l)p'{z) - Tp{z)) = J {t-^ -zt-{a + l)- rt'^) e^^e^^di 
= -J^td (e^We^*) - e^(*)e^*dt = te^^e^* 



due to the selection of the contour. This proves the first identity in (a). We leave the details 
of the proof of the remaining identities in (a) to the reader. 
Since 



V'jt) - V'js) _1112-a 

t - s t^s ^ f^s^ ^ ts^ ts 



1 



1 

t? 



the expression in (8.28) is obtained by direct substitution of those in (a) into the right hand 
side and comparison with (8.23). 

Furthermore, if F and S do not intersect, then 







d 



1 d 1 
+ 



r is L^* t — s ds t — s 



r 



d 



1 



Ir It, 



-V(s)-xs ^'it) + ^ V(t)+xt 



dsdt 



Ir I 



ds \t ~ s 

V(t)+xt ~^'i^) ~ ^ ^-V{s)-xs 
t - S 



^-V{s)-xs^^^^ 

dsdt 



where we have used integration by parts and (8.28). 

On the other hand, if the only intersection of F and S is at the origin, we can perform 
the same calculation with F and = S \ ?7e(0), where i7e(0) = {2; G C : \z\ < e}. Taking 
a posteriori e — t- and observing that all integrands are strongly vanishing at the origin, we 
arrive at (c) also in this case. 

Assume finally that there exists a point / such that F n S = {2:0} • Denote now 
Se = S\?7,(zo). We have 



■- lim 



= lim 



f [ e-^i.s)-xsTit)+^^Vit)+xt^^^^ 
^V(t)+xt -V'i^)-^ ^-Vis)-xs^^^^ 



r JE, 



r L 



]_^V{t)-V{s")+x{t-s") ^^V{t)-V{s')+x{t-s') 



dt 
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where s' and s" are the two points of intersection of E with the circle |z — zqI = e. We can 
deform the path of integration T in such a way that it forms a small loop around z = s' , 
picking up the reside of the integrand, and conclude that 



\p,q]{x) = 2m. 



□ 



With account of Lemma 8.3 we define g^, j = 1, 2, 3, as in (8.26), with paths Sj specified 
in Figure 15 and 



1 g«7ri 

^ 27ri ' ^ 27rz ' 

and conclude that condition (8.24) is satisfied. 

Let us turn to the equality (1.19); for that, let us define 



(8.30) 



K{x,y;z) = 



^V{t)-V{s)^xt-ys+{x-y) \og(z) 



ter Jses 



dtds 



(8.31) 



where the contours are as described in Theorem 1.2. A straightforward computation shows 
that 



d_ 



K{x,y;z) 



= {x- y)K{x, y; 1) = {x - y) 

z=\ JterJsei: 



^V(t)-V(s)^xt-ys 



dtds 
s-t' 



On the other hand, the change of variables t^t — log(2;) and s ^ s — log(2;) in the definition 
of K{x, y; z) yields 



K{x,y;z) 



^Vit-\ogiz))-V(s-\og(z))^xt-ys 



ter Jses 



dtds 



so that 



d_ 



K{x,y;z] 



=1 iter Jses * — ■s 
Thus, we get 

r I m^n£),v(.)-v,V'-..dM» = / 

JteT JseT, t — s 7(gr 7ses 



^V{t)-V{s)^xt-ys 



dtds 
s-t' 



Prom (8.28) where we plug in the previous identity, together with the value of C3 in (8.30) 
and the definition of p2 in (5.16), we obtain (1.19). 
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